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We review the B meson lifetime problems and nonspectator effects. The predictions of B 
meson lifetime ratios depend on four unknown hadronic parameters Si, B2, ei and C2, where 
B\ and B2 parametrize the matrix elements of color singlet-singlet four-quark operators and 
ei and C2 the matrix elements of color octet-octet operators. To understand contributions of 
the nonspectator effects to the B meson lifetime ratios, we derive the renormalization-group 
improved QCD sum rules for these parameters within the framework of heavy quark effective 
theory. 



I. INTRODUCTION 



A QCD-based operator-product-expansion (OPE) formulation for treatment of inclusive heavy hadron 
decays has been developed in past years [1]. The optical theorem tell us that the inclusive decay rates are 
related to the imaginary part of certain forward scattering amplitudes along the physical cut. Since, based 
, on the hypothesis of quark-hadron duality, the final state effects which are nonperturbative in nature, are 
C\| ' eliminated after adding up all of the states, the OPE approach thus can be employed for such the smeared or 
averaged physical quantities. In order to test the validity of (local) quark-hadron duality, it is very important 
to have a reliable estimate of the heavy hadron lifetimes within the OPE framework and compare them with 
ON ' experiment. 

. In the heavy quark limit, all bottom hadrons have the same lifetimes in the parton picture. With the advent 
Cl^ of heavy quark effective theory, which gives a systematic way in expansion of the initial heavy hadron, and the 
OPE approach for the analysis of inclusive weak decays, it is realized that the first nonperturbative correction 
to bottom hadron lifetimes starts at order l/m\. However, the 1/m^ corrections are small and essentially 
negligible in the lifetime ratios. The nonspectator effects such as W^-exchange and Pauli interference due 
to four-quark interactions are of order 1/mg, but their contributions can be potentially significant due to 
a phase-space enhancement by a factor of I671" 2 . As a result, the lifetime differences of heavy hadrons come 

I mainly from the above-mentioned nonspectator effects. 

II. DIFFICULTIES OF THE OPE APPROACH 

The world average lifetime ratios of bottom hadrons are [2] : 

r(B-) 



r(B° s ) 



= 1.07 ±0.03, 
0.94 ±0.04. 
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-^| = 0.79 ±0.05. (2.1) 

Since, to order l/m\, the OPE results for all of the above ratios are very close to unity [see Eq. (3.11) 
below] , the conflict between theory and experiment for this lifetime ratio is quite striking [3-6] . One possible 
reason for the discrepancy is that (local) quark-hadron duality may not work in the study of nonleptonic 
inclusive decay widths. Another possibility is that some hadronic matrix elements of four-quark operators 
are probably larger than what naively expected so that the nonspectator effects of order 167r 2 /m| may be 
large enough to explain the observed lifetime ratios. Therefore, one cannot conclude that (local) duality 
truly fails before a reliable calculation of the four-quark matrix elements is obtained [4] . 

Conventionally, the hadronic matrix elements of four-quark operators are evaluated using the factorization 
approximation for mesons and the quark model for baryons. However, as we shall see, nonfactorizablc effects 
absent in the factorization hypothesis can affect the B meson lifetime ratios significantly. To have a reliable 
estimate of the hadronic parameters B\, £?2, ei and 62 in the meson sector, to be introduced below, we will 
apply the QCD sum rule to calculate these unknown parameters. 



III. THEORETICAL REVIEW 

In this talk we will focus on the study of the four-quark matrix elements of the B meson. Before proceeding, 
let us briefly review the theory. Applying the optical theorem, the inclusive decay width of the hadron H b 
containing a b quark can be expressed as 

T(H b ^ X) = — ImJ d 4 x(H b \T{jC eS (x),jC eS (0)}\H b ) , (3.1) 
m Hb J 

where £ e fj is the relevant effective weak Lagrangian that contributes to the particular final state X. When 
the energy release in a b quark decay is sufficiently large, it is possible to express the nonlocal operator 
product in Eq. (3.1) as a series of local operators in powers of l/m b by using the OPE technique. In the 
OPE series, the only locally gauge invariant operator with dimension four, bipb, can be reduced to m b bb 
by using the equation of motion. Therefore, the first nonperturbative correction to the inclusive B hadron 
decay width starts at order l/mf } . As a result, the inclusive decay width of a hadron H b can be expressed 
as [1] 

G 2 F ml\V C Ku\ 2 1 / Xlv x {H b \b\g s a-Gb\H b ) 



- X) = ' ,- : r 5— [4 (H b \bb\H b) + d 
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+Ec x ( „) moiy b) +0{l/mt) , (3 _ 2) 



where a ■ G — ct^C*", Vqkm denotes some combination of the Cabibbo-Kobayashi-Maskawa parameters 
and cf reflect short-distance dynamics and phase-space corrections. The matrix elements in Eq. (3.2) can 
be systematically expanded in powers of l/m b in heavy quark effective theory (HQET), in which the 6-quark 
field is represented by a four- velocity-dependent field denoted by h^(x). To first order in l/m b , the 6-quark 
field b(x) in QCD and the HQET-field h£\x) are related via 



b(x) = e 



l + i 4 



2m b 



hW(x). (3.3) 



Applying this relation, one can replace b by the effective field in Eq. (3.2) to obtain 

(HMH b ) + G^ 

2m Hb 2m 2 b 2m\ W bJ 

{H b \b\g s a ■ Gb\H b ) 

= G Hb + 0{l/m b ) , (3.4) 
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where 



(H b \f& b) (iD ± )V\H b ) 



K Hi = - 



(H b \hi b) \g a a-Gh { v b) \H b ) 
2m H , 



(3.5) 



Note that here we adopt the convention D a = d a — ig s A a . The inclusive nonleptonic and semilcptonic decay 
rates of a bottom hadron to order 1 /m\ are given by [1] 



r NL (tf b ) = S^UlK 6 l 2 : 



1927T 3 



2cic 2 



2m Hb 

(al (x, 0, 0) + 0I o (x, x, 0)) (H b \bb\H b ) 
1 



mt 



(h (x, 0, 0) + h (x, x, 0)) (H b \bg s a ■ Gb\H b ) 



(3.6) 



where N c is the number of colors, the parameters a and (3 denote QCD radiative corrections to the processes 
b — > cud and 6 — + ccs, respectively [8], and 

sl( b) "T92^' cb| 



I o (x,0,0)(H b \bb\H b ) - — h{x,0,0)(H b \bg s a ■ Gb\H b ) 

771, 



(3.7) 



where r](x, xe, 0) with xi = (mt/mQ) 2 is the QCD radiative correction to the semileptonic decay rate and its 
general analytic expression is given in [7]. In Eqs. (3.6) and (3.7), /o,i,2 are phase-space factors (see e.g. [5] for 
their explicit expressions): 7i(x,0, 0) for b — > cud transition and li(x,x,0) for b — > ccs transition. Note that 
the CKM parameter V u d does not occur in TNh(Hb) and TsiX-Hb) when summing over the Cabibbo-allowed 
and Cabibbo-suppressed contributions. 

In Eq. (3.6) C\ and c 2 are the Wilson coefficients in the effective Hamiltonian 



<i/AB=l Gf 
H eS - 



V cb V: q (c 1 (»)OUri+c2(n)OUv)) 



V2 

+ v cb v: q ( Cl {n)oi{p) + C2 (m)o 2 c (^)) + 

where q — d,s, and 



+ h.c, 



2 =97m(1- 75)6 cf(l -7 5 )u. 



(3.6 



(3.9) 



The scale and scheme dependence of the Wilson coefficients ci^(^) are canceled out by the corresponding 
dependence in the matrix element of the four-quark operators Oi,2- That is, the four-quark operators in the 
effective theory have to be renormalized at the same scale \x and evaluated using the same renormalization 
scheme as that for the Wilson coefficients. 

Here we use the effective Wilson coefficients Cj which are scheme- independent [5] . 



ci = 1.149, 



c 2 = -0.325 . 



Using m b = 4.85 GeV, m c = 1.45 GeV, |V c6 | = 0.039, G B = 0.36 GeV 2 , G Ab = 0, K B « K Ab 
together with a = 1.063 and (3 — 1.32 to the next-to-leading order [8], we find numerically 



r(B-) 
r(B d ) 



r(B s ) 



r(A 6 ) 



= l + 0(l/mg), l±-lL = i + 0(l/ml), 3^ = 0.99 + 0(l/mg ) . 



r{B A ) 



r{B d ) 



(3.10) 
0.4 GeV 2 

(3.11) 
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It is evident that the 1/m 2 corrections are too small to explain the shorter lifetime of the A& relative to 
that of the Bd- To the order of l/m|, the nonspectator effects due to Pauli interference and W-exchange 
parametrized in terms of the hadronic parameters [4]: B\, B 2 , e\, e 2 , B, and r (see below), may contribute 
significantly to lifetime ratios due to a phase-space enhancement by a factor of 167r 2 . The four-quark 
operators relevant to inclusive nonleptonic B decays are 

O q v _ A = Ilj^l q L ^b L , 
O q s _ P = b R q L q L b R , 



T^_ A = b L ^t a q L q L ^t a b L , 
T g s _ P = b R t a q L q L t a b R , 



(3.12) 



where q Rt L = ^^-1- For the matrix elements of these four-quark operators between B hadron states, 
following, [4] we adopt the definitions: 

fk m B q 



1 

2m Bq 
1 

1 

1 

1 

2m Ab 
1 

2m A , 



(B q \0« v _ A \B g )=- 

■{B q \0%_ P \B q y. 

-{B q \T*_ A \B q )i 

-{B q \Tl_ p \B q y. 

-(A 6 |0^_ a |A 6 )e 
■(A b \T^_ A \A b ), 



fk m B q 



fk m B q 
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B, 

-B 2 

-£2 , 



48 

1,~ l x 1 

— (B+-) 

2 3 2mA b 



(A b \O q v _ A \A b ) 



(3.13) 



Under the factorization approximation, Bi = 1 and = 0, and under the valence quark approximation 
B = 1 [4]. 

The destructive Pauli interference in inclusive nonleptonic B~ decay and the W-exchange contributions 
to B° and B° s are [4] 1 



-r |Kd| 2 %s P cc(l - 2;) 2 < (1 + \x) [(-^-c 2 + 2c lC2 + JV c t|)Si + 2c 2 £l 



-(1 + 2s) (— c 2 + 2c lC2 + A c c 2 )5 2 + 2c 2 e 2 

J V r. 



-r |Kd| 2 r/„ S pccVl - 4z I (1 + ix) [(^c 2 + 2c lC2 + ^ c c 2 )B! + 2c 2 £l 



-(1 + 2s) {—c\ + 2c lC2 + N c c 2 2 )B 2 + 2c 2 e 2 



r^(B-) =r 77 nspcc (l-2;) 2 [(c?+^)(Bi+6ei)+6ciC2Bi] , 

r ann (s°) 



r |K s | 2 ?7n SP ocVl - 4a; j (1 + ^x) [(^c 2 + 2 Cl c 2 + N c c 2 2 )B x + 2c\e x 



1 The penguin- like nonspectator contributions to B a are considered in [10], but they are negligible compared to that 
from the current-current operators Oi and Oi introduced in Eq. (3.9). 
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-(1 + 2x) [(— c\ + 2c lC2 + N c c 2 2 )B 2 + 2c\e 2 



-r |K,| 2 % sp cc(l - xf |(1 + [(-^c? + 2c lC2 + iV c ^)Bi + 2c? ei 
-(1 + 2x) [(^-c? + 2cic 2 + Ar cC 2)B 2 + 2c 2 e 2 



(3.14) 



with 



ml 



where fs q is the B q meson decay constant defined by 

(0\qi^5b\B q {p)) = ij Bq Vp, ■ 
Likewise, the nonspectator effects in inclusive nonleptonic decays of the baryon are given by [4] 

r ann (A 6 ) = ir r/ nspcc r(l - xf (b(c\ + c\) - 2c lC2 ) , 

r!_ nt (A 6 ) = -ir [|K d | 2 (i - xf{\ + x) + \v cd \ 2 VT^] (b c \ - 2 Cl c 2 - jv cC £) . 

Using the values of c, in Eqs. (3.10), we obtain 

r ann (B d ) = T r/ nspcc (-0.0087B 1 + 0.0098B 2 - 2.28 £l + 2.58e 2 ) , 
T mt (B-) = T r ?nspec (-0.68 J B 1 + 7.10ei) , 

r ann (B s ) = T r/ nspcc (-0.0085 J B 1 + 0.0096^2 - 2.22ei + 2.50e 2 ) , 
T ann (A b ) = T ry nspcc r(0.59B + 0.31) , 
T int (A b ) = T Vnspec r{-0.30B - 0.097) . 

Therefore, to the order of l/«i|, the _B-hadron lifetime ratios are given by 



t(B~ 



r(B° d ) 
r(B° s ) 
r(B° d ) 



f. 



B 



185 MeV 



) 2 (o.043Bi + 0.0006^2 - 0.61ei + 0.17e 2 ) , 



= 1+ ( ttJth Tr V(~l-7x 10 _5 Bi + 1.9x 10" 5 B 2 - 0.0044ei + 0.0050 e 2 ), 
V185 MeV/ 



0.99 + 



( fB V 
V185 MeV / 



O.OOO6B1 + 0.0006B 2 



,185 MeV 
-0.15ei + 0.17e 2 - (0.014 + 0.0195)r 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



We see that the coefficients of the color singlet-singlet operators are one to two orders of magnitude smaller 
than those of the color octet-octet operators. This implies that even a small deviation from the factorization 
approximation ej = can have a sizable impact on the lifetime ratios. It was argued in [4] that the unknown 
nonfactorizable contributions render it impossible to make reliable estimates on the magnitude of the lifetime 
ratios and even the sign of corrections. That is, the theoretical prediction for t(B~) /r(B d ) is not necessarily 
larger than unity. In the next section we will apply the QCD sum rule method to estimate the aforementioned 
hadronic parameters, especially ej. 
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IV. THE QCD SUM RULE CALCULATION 



In HQET where the b quark is treated as a static quark, we can use the renormalization group equation 
to express them in terms of the operators renormalized at a scale Aqcd <C \i <C m&. Their renormalization- 
group evolution is determined by the "hybrid" anomalous dimensions [11] in HQET. The operators O v _ A and 
Ty_ A , and similarly O q s _ P and Tj_ p , mix under renormalization. In the leading logarithmic approximation, 
the renormalization-group equation of the operator pair (O, T) reads 



(4.1) 

where t — ^ ln(Q 2 / fi 2 ), Cp = (N 2 — 1)/2N C , and effects of penguin operators induced from evolution have 
been neglected. 

The solution to the evolution equation Eq. (4.1) has the form 

/ vr ^ W LT M 




T 




where 



Q D„, (4.2) 



= "-"I , (-) 



Lq = a s (u)/a s (Q) and (3o = -3- N c — | nj is the leading-order expression of the /3-function with nj being the 
number of light quark flavors. The subscript /x in Eq. (4.3) and in what follows denotes the renormalization 
point of the operators. Given the evolution equation (4.2) for the four-quark operators, we see that the 
hadronic parameters Bi and e t normalized at the scale are related to that at \i = 1 GeV by 

B i (m6)^1.54B i (/i)-0.41e i (/i), 

e^mb) ~ -0.090Si(/i) + 1.07 C< (/i) , (4.4) 



with fi = 1 GeV, where uses have been made of a s {m z ) = 0.118, = 333 MeV, nib = 4.85 GeV, 
m c = 1.45 GeV. The above results (4.4) indicate that renormalization effects are quite significant. 

It is easily seen from Eqs. (4.2) and (4.3) that the normalized operator D\ (or D2) is simply multiplied 
by L/q 1 " 2 ^ (or 1) when it evolves from a renormalization point /1 to another point Q. In what follows 2 , wc 
will apply this property to derive the renormalization-group improved QCD sum rules for Dj at the typical 
scale /1 = 1 GeV. We define the new four-quark matrix elements as follows 

1 { B q \D®M\B 9 ) = ^^&M, (4.5) 



2m B , 

where the superscript (i) denotes (V — A) four-quark operators for i = 1 and (S — P) operators for i = 2 



and dj satisfy 




Bi — j€i 



(4.6) 



2 In the sum rule calculation, the factorization scale fj, cannot be chosen too small, otherwise the strong coupling 
constant a a would be so large that Wilson coefficients cannot be perturbatively calculated. 
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Since the terms linear in four-quark matrix elements are already of order 1/mf, we only need the 
relation between the full QCD field b(x) and the HQET field hi b \x) to the zeroth order in l/m b : 
b(x) = e - lm " v - x {h { y ] (x) + 0(l/m b )}. In the following, within the framework of HQET, we apply the 
method of QCD sum rules to obtain the value of the matrix elements of four-quark operators. We consider 
the three-point correlation function 

ifi'Vo/) = * 2 / dxdye^~^-y(0\T{[q(x)T a hi h \x)} Df\o) [q(y)T p h^(y)^}\0) , (4.7) 



ce,/9 

where the operator is defined in Eq. (4.3) but with b — > and T a is chosen to be v a "f 5 (some further 
discussions can be found in [12]). 

The correlation function can be written in the double dispersion relation form 

, a f f ds ds' D «(0 
IV, M =J J—y^fP' • (4-8) 

The results of the QCD sum rules are obtained in the following way. On the phenomenological side, which 
is the sum of the relevant hadron states, this correlation function can be written as 

where A is the binding energy of the heavy meson in the heavy quark limit and ellipses denote resonance 
contributions. The heavy-flavor-independent decay constant F defined in the heavy quark limit is given by 

(0|g7 M 7 5 4 b) \B(v)) = iF(p)v^ . (4.10) 

The decay constant F(ji) depends on the scale /z at which the effective current operator is renormalized and 
it is related to the scale-independent decay constant fs of the B meson by 

F{m b ) = f B y/mE. (4.11) 

On the theoretical side, the correlation function can be alternatively calculated in terms of quarks and 
gluons using the standard OPE technique. Then we equate the results on the phenomenological side with 
that on the theoretical side. However, since we are only interested in the properties of the ground state at 
hand, e.g., the B meson, we shall assume that contributions from excited states (on the phenomenological 
side) are approximated by the spectral density on the theoretical side of the sum rule, which starts from 
some thresholds (say, Uij in this study). To further improve the final result under consideration, we apply 
the Borel transform to both external variables u and u'. After the Borcl transform [13], 

B[nJg%V)]= lim lim J-(-u/) m+1 [ / 7 ] m (-w) n+1 [^-] n nfi < ' ' (w, w') , (4.12) 

LO— ^OO 

r fixed 



-^jr fixed 

the sum rule gives 



F 2 (m b )^( M ) e _ A/tle _ A/t2 ^) = n-* ds rfs/e _ (s/tl+s7t2)/9QC D ; (4 13) 

16 Jo Jo 

where Wij is the threshold of the excited states and pQ CD is the spectral density on the theoretical side of 
the sum rule. Because the Borel windows are symmetric in variables t\ and t2, it is natural to choose t\ =t<i- 
However, unlike the case of the normalization of the Isgur-Wise function at zero recoil, where the Borel mass 
is approximately twice as large as that in the corresponding two-point sum rule [14], in the present case of the 
three-point sum rule at hand, we find that the working Borel windows can be chosen as the same as that in 
the two-point sum rule since in our analysis the output results depend weakly on the Borel mass. Therefore, 
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we choose t\ = ti = t. By the renormalization group technique, the logarithmic dependence a s ln(2i//z) can 
be summed over to produce a factor like [a s (/x)/a s (2i)] 7 . After some manipulation we obtain the sum rule 
results: 



F 2 (m b )F 2 ( M ) 2A/t 
16 



q,(2f) 
a s (/i) 



7^ 



where 



OPE e 



1 - 2(5 



a a (2t) ' 



1 - 2(5 



OPE i; , , ^ i(OPE)i pt;i>j , 



L -9/(2 /3o) Q 



1 

16 



OPE Biil - f OPE £M 
i OPE s . 2 + OPE ei)2 



fea • Gq)_ t{l e _ wij/t) + i^) t 2 (1 _ e -. 1 ,/*)2 



8^ 2 



0PE £2 ~0(« s ), 



(4.14) 



(4.15) 



with 



(OPE) 2 pt;i J j 



a s /17 47r' „, s\ 

1 + v(y + -- 21n i) 



\ 1 + —) {qq)+ }■ 



(4.16) 



For reason of consistency, in the following numerical analysis we will neglect the finite part of radiative one 
loop corrections in OPE^ . and OPE Ci . (and in Eq. (4.19)). The parameter S in (4.14) is some combination 
of the (3 functions and anomalous dimensions (see Eq. (4.2) of [15]) and is numerically equal to —0.23. The 
relevant parameters normalized at the scale t are related to those at [x by [15,13] 



F{2t)=F{^ 

(qq)2t = (qq)^ 
(g s qa ■ Gq) 2t = 



a s (2i)\-2//3o 1 



1 j "*( 2 *) ' 



/ q a (2t) \-Vft> 
\a s (n)) 

{9sqa ■ Gq) » ■ {^M ) 



2/(3/3o ) 



(a s G 2 } 2t -- 

where (• • •) stands for (0| ■ ■ ■ |0) and [13] 



(a s G 2 ) 



(4.17) 



(qq)v=i Gov = -(240 McV) 3 , 

(« S G 2 ) M=1 GoV = 0.0377 GeV 4 , 

{qgsc^G^q)^! Gov = (0.8 GeV 2 ) x (qq)^ GcV . 



(4.18) 



Some remarks are in order. First, in Eqs. (4.14) and (4.15). OPEs i is obtained by substituting D^ 1 by 
O v and it can be approximately factorized as the product of (OPE)2 P t;i,j with itself, which is the same as the 
theoretical part in the two-point F(/j,) sum rule [14,15]. In the series of (OPE)2 P t ; ij, we have neglected the 
contribution proportional to (qq) 2 ■ (More precisely, it is equal to a s (g<7) 2 7r/324; see Ref. [14].) Nevertheless, 
the result of (OPE)^ in Eq. (4.15) is reliable up to dimension six, as the contributions from the (qq) 2 terms 
in (OPE) 2 pt ; ij are much smaller than the term (1 + a s / 'tt) 2 {qq) 2 / '16 that we have kept [see Eq. (4.16)]. 
Second, in (OPE)^ the contribution involving the gluon condensate is proportional to the light quark 
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mass and hence can be neglected. Third, OPE ei is the theoretical side of the sum rule, and it is obtained 
by substituting D-f l) by T v . Here we have neglected the dimension-6 four-quark condensate of the type 
(qTX a q qTX a q). Its contribution is much less than that from dimension- five or dimension- four condensates 
and hence unimportant (see [16] for similar discussions). It should be emphasized that nonfactorizablc 
contributions to the parameters Bi arise mainly from the O v — T v operator mixing. 

In the following, we compare our analysis with the similar QCD sum rule studies in [16] and [6]. First, 
Chernyak [16] used the chiral interpolating current for the B meson, so that all light quark fields in his 
correlators are purely left-handed. As a result, there are no quark-gluon mixed condensates as these require 
the presence of both left- and right-handed light quark fields. Instead, the gluon condensate contribution 
enters into the e± sum rule with an additional factor of 4 in comparison with ours; thus their OPE ei is in 
rough agreement with ours. Second, our results for OPE £i are very different from that obtained by Baek et 
al. [6]. The reason is that their results are mixed with the 1 + to 1 + transitions. Also a subtraction of the 
contribution from excited states is not carried out in [6] for the three-point correlation function, though it 
is justified to do so for two-point correlation functions. Indeed, in the following analysis, one will find that 
after subtracting the contribution from excited states, the contributions of OPE ei are largely suppressed. 
Furthermore, as in the study of the B meson decay constant [14], we find that the renormalization-group 
effects are very important in the sum rule analysis. Moreover, e, at /i = are largely enhanced by 
renormalization-group effects. 

The value of F in Eq. (4.14) can be substituted by 



F 2 (n)e 



-A/t _ 



~a s {2t)~ 


4 




u 


_ a s (fi) _ 




1 28 aA}l) 

IT 





3 



1 + 



a s (2t)/17 4tt 

e ' — i -t- 

2a s (2t)\._ (qg s (T-Gq) 



TT \ 3 9 t) 



)(«)»+ y , (4-19) 



which is from the two-point sum rule approach [15]. Next, to determine the thresholds Wjj we employ the 
B meson decay constant fs = (185 ± 25 ± 17) MeV obtained from a recent lattice-QCD calculation [17] and 
the relation [18] 

f B = *W(l- 2 -^) ( 1 - ( °- 8 ~ L1) GCV ) , (4.20) 

that takes into account QCD and 1/mfe corrections. Using the relation between F{m^) and F(fi) given by 
Eq. (4.17) and m b = (4.85 ± 0.25) GeV, we obtain 

F(n = 1 GeV) S (0.34 ~ 0.52) GcV 3/2 . (4.21) 

Since the A parameter in Eq. (4.19) can be replaced by the A sum rule obtained by applying the differential 
operator t 2 d\n jdt to both sides of Eq. (4.19), the F(p) sum rule can be rewritten as 

d 

F 2 {ijl) = (right hand side of Eq. (4.19)) x exp[t — ln(right hand side of Eq. (4.19))] , (4.22) 

which is A-free. Then using the result (4.21) as input, the threshold uj in the F(/j,) sum rule, Eq. (4.22), 
is determined. The result for ojq is 1.25 — 1.65 GeV. A larger F(p = 1 GeV) corresponds to a larger uq. 
The working Borel window lies in the region 0.6 GeV < t < 1 GeV, which turns out to be a reasonable 
choice. Substituting the value of ojo back into the A sum rule, we obtain A = 0.48 — 0.76 GeV in the Borel 
window 0.6 GeV < t < 1 GeV. This result is consistent with the choice mj = (4.85±0.25) GeV, recalling 
that in the heavy quark limit, A = tub — To extract the d v ^ sum rules, one can take the ratio of 
Eq. (4.19) and Eq. (4.14) to eliminate the contribution of _F 2 /exp(A/i). This means one has chosen the 
same A both in Eq. (4.19) and Eq. (4.14). Since quark-hadron duality is the basic assumption in the QCD 
sum rule approach, we expect that the same result of A also can be obtained using the A sum rules derived 
from Eq. (4.14) (see [13] for a further discussion). This property can help us to determine consistently the 
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threshold in 3-point sum rule, Eq. (4.14). Therefore, we can apply the differential operator i 2 <91n /dt to both 
sides of Eq. (4.14), the d v ^ sum rule, to obtain new A sum rules. The requirement of producing a reasonable 
value for A, say 0.48 — 0.76 GeV, provides severe constraints on the choices of Wij. With a careful study, we 
find that the best choice in our analysis is 

= -0.02 GcV + ^o, wi,2 = -0.5 GeV + ui , u 2 , 2 = -0.22 GeV + lo • (4.23) 

Applying the above relations with oj = (1-25 <~ 1.65) GeV and substituting F(fj,) in Eq. (4.14) by (4.19), 
we study numerically the d v ^ sum rules. In Fig. 1, wc plot and e\ as a function t, where B\ = 

8< W /9 + 2d v 2 (l) /3, and t\ = -4^ (i) /27 + 8d v 2 (l) /9. The dashed and solid curves stand for B? and e V, 
respectively, where we have used uiq = 1.4 GeV (the corresponding decay constant is Jb = 175 ~ 195 MeV 
or F(pb) = 0.405 ± 0.005 GeV 3 ' 2 ). The final results for the hadronic parameters Bi and ej are (see Fig. 2) 

BKfi = 1 GeV) = 0.60 ± 0.02, B^(ji = 1 GeV) = 0.61 ± 0.01, 

el(fj,= l GeV) = -0.08 ± 0.01, e^/U = 1 GeV) = -0.024 ± 0.006. (4.24) 

The numerical errors come mainly from the uncertainty of luq = 1.25 ~ 1.65 GeV. Some intrinsic errors of 
the sum rule approach, say quark-hadron duality or a s corrections, will not be considered here. 
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FIG. 1. B?(M) and e^M) as a function t, where E° = U\ w /9 + 2(^ w /3, and = -4d? w /27 + 8d2 tl 79. The 
dashed and solid curves stand for and e", respectively. Here we have used u>o = 1.2 GeV and Eq. (4.23). 



Substituting the above results into Eq. (4.4) yields 

0.96 ± 0.04 + 0(l/m b ), 
-0.14 ± 0.01 + 0(l/m b ), 

It follows from Eq. (3.19) that 

r(B-) 



Bi(m b ) 
ei{m b ) = 



r(B d ) 
r(B s ) 
r(B d ) 

r(B d ) 



1.11 ±0.02, 



1 



0.99 



(l^cv)^°- 007 + - 020i?)r ' 



B 2 (m b ) = 0.95 ± 0.02 + 0(l/m b ) , 
£2(m 6 ) = -0.08 ± 0.01 + 0(l/m b ) . 



(4.25) 



(4.26) 



to the order of 1/ jrifo. Note that we have neglected the corrections of SU(3) symmetry breaking to the 
nonspectator effects in T(B s )/r(Bd). We see that the prediction for T(B~)/r{Bd) is in agreement with 
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the current world average: t(B~)/t(B ( i) =1.07± 0.03 [2], whereas the heavy-quark-expansion-based result 
for t(B s )/t(Bci) deviates somewhat from the central value of the world average: 0.94 ± 0.04. Thus it is 
urgent to carry out more precise measurements of the B s lifetime. Using the existing sum rule estimate 
for the parameter r [20] together with B = 1 gives T(Af,)/r(J5 t ;) > 0.98. Therefore, the 1/mf nonspectator 
corrections are not responsible for the observed lifetime difference between the A& and B4. 



V. CONCLUSIONS 



The nonspectator effects can be parametrized in terms of the hadronic parameters Bi, B 2l ei and e 2 [21], 
where B\ and B2 characterize the matrix elements of color singlet-singlet four-quark operators and e\ and 
€2 the matrix elements of color octet-octet operators. In OPE language, the prediction of B meson lifetime 
ratios depends on the nonspectator effects of order \6ir 2 /m\ in the heavy quark expansion. Obviously, the 
shorter lifetime of the A b relative to that of the Bd meson and/or the lifetime ratio t(B s )/t(B,i) cannot be 
explained by the theory so far. It is very likely that local quark-hadron duality is violated in nonlcptonic 
decays. 

As emphasized in [5], one should not be contented with the agreement between theory and experiment 
for the lifetime ratio T{B~)/r{Bd). In order to test the OPE approach for inclusive nonleptonic decay, it is 
even more important to calculate the absolute decay widths of the B mesons and compare them with the 
data. From (3.6), (3.7), (4.25) and considering the contributions of the nonspectator effects, we obtain 

T tot (B d ) = (3.61™) x 10- 13 GeV, 

rtot(B-) - (3.34™) x 10- 13 GeV, (5.1) 

noting that the next-to-leading QCD radiative correction to the inclusive decay width has been included. The 
absolute decay widths strongly depend on the value of the b quark mass. The problem with the absolute decay 
width T(B) is intimately related to the B meson semileptonic branching ratio Bsl- Unlike the semileptonic 
decays, the heavy quark expansion in inclusive nonleptonic decay is a priori not justified due to the absence 
of an analytic continuation into the complex plane and hence local duality has to be invoked in order to 
apply the OPE directly in the physical region. 

To conclude, we have derived in heavy quark effective theory the renormalization-group improved sum rules 
for the hadronic parameters B\, B 2 , ei, and e 2 appearing in the matrix element of four-quark operators. The 
results are Bi(m 6 ) = 0.96 ±0.04, B 2 (m b ) = 0.95 ±0.02, ei(m 6 ) = -0.14 ±0.01 and e 2 (m b ) = -0.08 ±0.01 
to the zeroth order in 1/m;,. The resultant B-meson lifetime ratios are t(B~)/t(B4) = 1.11 ± 0.02 and 
T{B s )/r{B d ) w 1. 
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